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<D 

^ We compute the total cross section for the production of a light CP even Higgs 

boson within the framework of supersymmetric QCD up to next-to-next-to-leading 

i-H order. Technical subtleties in connection to the evaluation of three-loop Feynman 

integrals with many mass scales are discussed in detail and explicit results for the 

qq counterterms of the evanescent couplings are provided. The new results are applied 

I/") to several phenomenological scenarios which are in accordance with the recent dis- 

covery at the LHC. In a large part of the still allowed parameter space the K factor 
of the supersymmetric theory is close to the one of the Standard Model. However, 

^vq for the case where one of the top squarks is light, a deviation of more than 5% in the 

I next-to-next-to-leading order prediction of the cross section can be observed where 

at the same time the predicted Higgs boson mass has a value of about 125 GeV. 

PACS numbers: 12.38.Bx 12.60.Jv 14.80.Da 



1 Introduction 



Recently the experiments ATLAS and CMS at the Large Hadron Collider (LHC) at 
CERN reported on the discovery of a new particle with a mass of about 125 GeV [T],[2J . 
Although not all of its properties have been determined yet the new particle seems to be 
in accordance with a Higgs boson, maybe the one of the Standard Model (SM), maybe 
one corresponding to extensions of the SM. In this paper we consider a supersymmetric 
extension of the SM and report about next-to-next-to-leading order (NNLO) corrections 
to the production cross section of the lightest CP even Higgs boson in this model. One 
of the aims is to quantify the size of the supersymmetric corrections. First results in this 
direction have been reported in Ref. [31. We will significantly extend these considerations 
and furthermore provide details on the calculation. 

The radiative corrections both for inclusive and exclusive processes are summarized in the 
reports (2) and (H], respectively. The numerical results for the SM Higgs boson production 
cross section provided in these papers are based on calculations to the gluon fusion process 
which have been performed up to next-to-leading order (NLO) [6h9| for arbitrary values of 



the top quark and Higgs boson mass and to NNLO in Refs. [10 -12 adopting the infinite 
top quark mass approximation. NNLO results for finite top quark masses have been 



obtained in Ref. 13 - 18 



Within the SM the bottom quark contributions are small and NLO terms [9] are sufficient 
considering the size of the corrections induced by the top quark Yukawa coupling. This 
is different in the case of supersymmetry (SUSY) where the factor m& is accompanied 
by tan (3, the ratio of the vacuum expectation values of the two Higgs doublets, and can 
thus lead to large contributions to the cross section for the production of the lighter 



CP even Higgs boson. Explicit analytical NLO results can be found in Refs. [19 20 



(see also Ref. 21 for numerical results). The contribution from the top quark/top squark 
sector is available already since a few years. Within the effective-theory framework, which 
constitutes an excellent approximation for Higgs boson masses below about 200 GeV, the 



production cross section has first been computed in Refs. [22 23 and has been confirmed 
analytically [24] and numerically [2l] (see also Ref. [25])- In Refs. 26,27 the squark loop 
contributions to Higgs boson production in the Minimal Supersymmetric Standard Model 
(MSSM) have been computed without assuming any mass hierarchy. 

The first supersymmetric QCD corrections at NNLO have been computed in Ref. [3] (see 



also Refs. 23 , 28 for estimates prior to this work) for a degenerate mass spectrum. In 
this work these considerations are extended to allow for a wide range of phenomenologi- 
cally interesting scenarios. The three-loop results for the matching coefficients which are 



presented in this paper have been confirmed by an independent calculation in Ref. 29 
Note that the calculations in the paper are performed in the framework of the MSSM, 
however, they are in principle not limited to the minimal version and could, e.g., also 
be applied to the so-called Next-to-Minimal Supersymmetric Standard Model (NMSSM) 
where an additional singlet Higgs field is added as compared to the MSSM. 
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In Ref. [3] all supersymmetric masses have been identified. Furthermore, some of the 
technical difficulties in connection to the interplay of dimensional reduction (DRED) 
and supersymmetry have been addressed in this reference. In this context it has been 
very useful to evaluate the QCD corrections within DRED and check that the proper 
renormalization of the evanescent couplings leads to the known analytical expressions after 
re-expressing the strong DRED coupling by its counterpart in dimensional regularization 
(DREG). In this paper we build on the developed formalism and evaluate the Higgs-gluon 
coupling for various hierarchies which are of phenomenological relevance. In particular, 
we provide the counterterms of the evanescent couplings which could be relevant also for 
other calculations. 

The outline of the paper is as follows: in the next section we briefly review the formalism 
and discuss the subtleties at NLO (where an exact calculation has been performed) and 
NNLO. Moreover the quality of the approximations applied at NNLO are tested at the 
NLO level. We discuss at length the renormalization procedure. Numerical results for 
the Higgs-gluon coupling are also presented in this Section. Afterwards we compute in 
Section [3] predictions for the total Higgs boson production cross section and consider 
different scenarios: m™ ax which maximizes the prediction of the lighter CP even MSSM 
Higgs boson and a scenario where one of the top squarks is light whereas the reminder of 
the mass spectrum is heavy. Our conclusions are presented in Section |4| 



2 Higgs-gluon coupling in the MSSM 
2.1 Effective theory 

The theoretical framework needed for the computation of NNLO corrections to the pro- 
duction of a light CP even Higgs boson has already been described in Ref. (3j. For 
completeness we repeat the most important issues and formulae. Starting point is the 
effective five-flavour Lagrange density 

£y,efr = -^-C°C? + £qcd > W 

where Cq CD is the usual QCD part with five massless quarks. 0° denotes the light CP- 
even Higgs boson of the MSSM, v° is given by a/ (^i) 2 + ("W^) 2 where v\ and v% are the 
vacuum expectation values of the two Higgs doublets, and is the gluonic field strength 

tensor constructed from fields and couplings already present in c!q CD . The superscript 
"0" denotes bare quantities. Note that the renormalization of <p° /v is of higher order in 
the electromagnetic coupling constant. In Eq. Q the operator is given by 

Ol = \G%G ^. (2) 
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Figure 1: Sample diagrams contributing to G\ at one, two and three loops originating 
from the coupling of the Higgs boson to top quarks and top squarks. The symbols t, U, 
g, g, h and e denote top quarks, top squarks, gluons, gluinos, Higgs bosons and e scalars, 
respectively. 

It contains the coupling of the Higgs boson to two, three and four gluons. Heavy degrees 
of freedom only contribute to the coefficient function C°. In the SM this concerns only 
the top quark whereas for supersymmetric QCD also squarks and gluinos contribute. 

From the technical point of view C° is computed from vertex diagrams involving two 
gluons and a Higgs boson. In the practical calculation it is convenient to subdivide the 
contributing Feynman diagrams into two classes. The first one contains all contributions 
where the external Higgs boson couples to a top quark or top squark (see Fig. [T] for sample 
diagrams) whereas all diagrams with a coupling of the Higgs boson to a super partner 
of the five lighter quark flavours constitute the second class (see Fig. [2j. Note that the 
diagrams of the latter class lead to corrections to C° which are suppressed by M§/m~ 
whereas the Higgs-top-squark coupling generates contributions proportional to 1, to|/to~, 
m 4 /i susy /m? and M§/m| (where fi susy is the Higgs-Higgsino bilinear coupling from the 
super potential). Although the relative suppression is formally only M\jvn 2 t all Mf/m? 
terms are numerically much smaller than the leading contribution. In our calculation we 
consider all squarks except the top squark to be degenerate and we furthermore assume 
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Figure 2: Sample diagrams contributing to G\ at one, two and three loops originating from 
the coupling of the Higgs boson to the squarks corresponding to the light quarks. The 
symbols q, cji, g, g, c, h and e denote light quarks, the corresponding scalar super partners, 
gluons, gluinos, ghosts, Higgs bosons and e scalars, respectively. a q is an auxiliary particle 
used for the implementation of the four-squark vertex. 



that the corresponding quarks are massless. Thus the squark contributions from the 
first two families cancel and only the contribution from bottom squarks remains. In the 
following we nevertheless use the notion "squark contribution" for diagrams from Fig. [2j 

In this paper we compute all contributions from the top quark/top squark and light 
quark/squark sector in the limit of vanishing light quark masses up to three loops. For 



the corresponding Feynman rules we refer to the Appendices of Refs. 23 and 20 . Note 



that the couplings to the super partners of the light quarks has been implemented by 
averaging over the up, down, strange, charm and bottom squark contribution. Since for 
large values of tan j3 bottom quark diagrams might give large contributions they will be 



included in our numerical prediction up to NLO 19 20 as we will describe in Section 3.1 



Since C° constitutes the effective Higgs-gluon coupling it is sufficient to evaluate the 
corresponding Feynman integrals in the limit of vanishing external momenta. This is 
conveniently achieved by applying a projector (see, e.g., Eqs. (11) and (12) in Ref. |3l) on 
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one of the two tensor structures and by performing a naive Taylor expansion afterwards. 
As a result one ends up with vacuum integrals involving several mass scales in the case 
of supersymmetry. Up to two-loop order such integrals can be computed exactly, the 
corresponding algorithm and necessary formulae have been provided in Ref. [30] almost 
20 years ago. Since an analogue algorithm at three-loop order is not available to date one 
has to rely on special cases combined with expansions in order to construct approximations 



to the exact formulae. A general outline of such an approach has been discussed in Ref. 31 



where three-loop corrections to the lightest Higgs boson mass have been computed. In this 



paper we refine the approximation method which is described in Subsection 2.3 However, 
we first discuss in the next subsection the renormalization procedure. 



2.2 Renormalization and evanescent couplings 



Besides the technical difficulties, which are described in Subsection 2.3 also several field 



theoretic challenges have to be solved. They are mainly connected to the fact that for the 



calculation of C° it is convenient to use DRED 32 35 whereas the calculations within 
the effective theory of Eq. Q are performed with the help of DREG. The LO calculation 
is not affected by this issue. Also at NLO it is possible to proceed in a naive way. At 
NNLO, however, it is necessary to properly renormalize all occurring parameters and to 
use the relations between their DRED and DREG definition. In Ref. [3] the procedure 
has been described for a degenerate SUSY spectrum, the modifications necessary to cover 
the most general case are straightforward and briefly touched in the following. 

We renormalize all parameters in the DR scheme except M e and the coupling of the Higgs 
boson to e scalars. M e is renormalized on-shell with the condition M e = 0. Note that 
we take this limit only in the final result since in intermediate steps 1/M e terms occur 
which have to be canceled by proper counterterms. The two-loop DR counterterms for a s 



and m t are independent of all occurring scales and can be found in Refs. 31 , 36] (see also 
references therein). Due to the quadratic divergence the renormalization constants for 
the squark masses and the mixing angle 9 t , which we also need up to two loops, depend 
on mass ratios which is known exactly, even up to three loops [36] . The gluino mass 
counterterm is only needed to one-loop order. As far as wave function renormalization of 
external particles is concerned only the one for the external gluons contributes. Actually, 
for our calculation only those contributions to the gluon two-point function have to be 
considered which contain heavy particles. This is conveniently formulated in terms of a 



bare decoupling constant of the gluon field 37 which is usually denoted by £f. Note that in 



our case (® has to accomplish both the decoupling of all supersymmetric particles and the 
top quark and the transition from DRED to DREG. For the degenerate mass case results 
can be found in Ref. |3|. The analytical expressions needed for the calculations performed 



in this paper can be found in electronic form in 38 . Note that C° is multiplicatively 
renormalized with I/C3. 

A non-standard part of the renormalization procedure concerns the so-called evanescent 
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couplings which emerge through radiative corrections. In our case there are two operators 
which define the following Lagrange density (see also discussion in Ref. [2~T]) 



0\ 2 _0,o_0,a 



(3) 



where e 0,a denotes the bare e scalar field and the dimensionful quantity A° mediates the 
coupling of the Higgs boson to e scalars. The on-shell e scalar mass counterterm is needed 
to one- loop accuracy and is given by (see, e.g., Ref. |39|) 



o\2 



;m £ °) 



5M 




rnj 



3 2 
2™* 



m 7 



+ 



+ — 

8 2 



si: 



6 8 8 2 2 
+ O (M £ 4 ) + O (e 2 ) , (4) 

with l x = ln(/x^/ 'm 2 ,) and l e = \a.{^ 2 R /Mf) where /ir is the renormalization scale, ni — 5 
counts the number of massless quarks. Note that for our calculation the terms of O(M^) 
have to be kept until the inverse powers in the e scalar mass are canceled (cf. discussion 
below) . 

The counterterm for A £ is defined through 



5Kt + A s 



(5) 



and is needed to two-loop accuracy. Sample Feynman diagrams contributing to 5A £ 
are shown in Fig. [3j We apply the renormalization condition A £ = 0. The analytical 
expressions corresponding to the three hierarchies are quite lengthy — in particular for 
those cases where 8t ^ — and can be found on the webpage 38|. To illustrate the 

^ 2 



structure of the result we present in the following 5K £ for the hierarchy (hi) which is 



defined in Eq. (14). It is convenient to split the result into a contribution from the top 



quark/top squark sector and from the light quark/squark sector. Thus the leading terms 
read 




-2M 2 z s a 
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Figure 3: One- and two-loop sample diagrams contributing to 5A £ . The notation is 
adopted from Figs. [I] and [2} 
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where m susy = = = = m^, c x = cos(x), s x = sin(x) and t x = tan(:r). Here a is 
the mixing angle between the weak and the mass eigenstates of the neutral scalar Higgs 
bosons and tan is the ratio of the vacuum expectation values of the two Higgs doublets. 

A closer look to Eq. ^ shows that the terms proportional to M| exactly cancel in the sum 

of (sA^pj and (sA^j . Thus for a massless top quark there is no need to introduce a 
counterterm for A £ since the corresponding contributions cancel within one family. Note, 
however, that the M§ terms are needed to obtain separately finite results for the top and 
light quark sector. 

In the supersymmetric limit the appearance of the evanescent coupling 5X £ is forbidden by 
supersymmetry and thus there will be no radiative corrections to a vanishing hee coupling. 
Both 5A £ i and 5A e ^ vanish independently. Any non vanishing radiative corrections to the 
hee coupling is generated by the presence of soft SUSY breaking terms in the Lagrangian. 

It is interesting to note that without introducing 5A £ the terms in C\ proportional to M§ 
are not finite in the limit M £ — » 0. Rather they behave like M|/M 2 at two- loop order and 
even like M|/M^ at three loops which is due to diagrams like the third one in the third 
row of Fig. [TJ The inverse powers in M £ only cancel after the contribution from SA £ is 
added and thus the limit for vanishing e scalar mass can be taken in the final expression. 
The cancellation works independently for the top (Fig. [T| and light quark sector (Fig. [2|. 
In Eq. (|6j) we also list the contributions of order M 2 /m 2 and M 2 /m 2 usy which, however, 
do not contribute to the final result. In intermediate steps also m 2 /M 2 terms are present 
which, however, cancel among the diagrams involving top quark and top squark loops. 

From this discussion it becomes clear that the 1/M" terms could be avoided if a matching 
to an effective Lagrangian defined within DRED would be performed. However, in that 
case it is likely that more operators have to be considered in the effective Lagrangian 
which have not yet been studied in the literature. Alternatively the calculations could 
probably also be performed in the limit where M £ is the largest scale. Since the final result 
has to be independent of M £ one could perform the limit M £ — > oo. As a consequence, 
the inverse powers in M £ drop out by construction and the positive powers have to cancel 
in the final result, in analogy to the inverse powers in our approach. 

The next step concerns the transformation of the DR strong coupling in the full theory 
(ats ) to the MS version in five-flavour QCD («i 5) ). This is again achieved with the 
help of the corresponding decoupling constant defined through 

«i 5) = (a s af^K (7) 
Results for ( as can again be found in Refs. (3) for the degenerate mass case. Exact results 



can be obtained from the calculations performed in Ref. 40 but also from 29 where an 
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renormalization constant 


explanation 


loop order 


Z as 


strong coupling constant 
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H 


top squark mass, i = 1, 2 


2 


V 

Ttlq 


squark mass 


2 


y 


gluino mass 


f 


ml 


e scalar mass 


1 




hee coupling 


2 


Z n 


operator renormalization 


2 


decoupling constant 


explanation 


loop order 


AO 

S3 


decouples heavy particles form gluon field 


2 


Ca s 


relates af 1 to ai 3 * 300 -* 


2 



Table 1: Counterterms and decoupling constants entering the three-loop calculation of 



explicit result is provided. Note that it is important to have at hand the one-loop term 
of ( as up to order e since at this step still contains poles originating from the light 
degrees of freedom. They are removed with the help of the renormalization constant 
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7T e 





2 



/ ai 5 '' \ 1/51 19n 



via 



6/ + - e rirl' ^ 



^ = T^C'i , (9) 



where /3^ 5 ) = -(a { s ] /Tif(3f ] + ... is the five-flavour QCD beta function with /3^ l) = 
11/4 — ni/6 and in the second line the perturbative expansion up to two loops is displayed. 

For convenience we list in Tab. [T] all counterterms and decoupling constants needed for 
our three-loop calculation together with the required loop order. 

Let us remark that at LO and NLO the contributions to C\ from diagrams as shown in 
Fig. [T] and in Fig. [2] are separately finite after taking into account all relevant countert- 
erms. At NNLO, however, the contribution from the top sector contains still 1/e poles 
proportional to M§ / m^ usy (where m susy denotes a squark mass) which only cancel in com- 
bination with the result where the Higgs boson couples to a super partner of one of the 
five light quarks (cf. Fig. [2]). Sample diagrams responsible for such poles are given in 
Fig. [I] (first diagram in second row) and Fig. [2] (third diagram in second row). 
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At this point one has arrived at a finite result for C\ expressed in terms of a 
its perturbative expansion in the form 

2 



(5) 
s • 



Ci 



(5) 



.(0) 



Oft) + 



a 



(5) 



7T 



Oft) + 




We write 



(10) 



where we have introduced the renormalization scale fih ("h" stands for "hard"). The 
one-loop expression is collected in 111 the SM and the MSSM it is given by 
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mtHsusy — 2 2~ ' ( 12 ) 

V m| mi 

respectively, where c x = cos(x), s x = sin(x), t x = tan(x), sw is the sine of the weak 
mixing angle 9w, /^susy is the Higgs-Higgsino bilinear coupling from the super potential 
and tan is the ratio of the vacuum expectation values of the two Higgs doublets. The 
term proportional to Mf/m~ arises from the one-loop diagrams involving squarks, all 



other terms origin from the top quark and top squark loops. Note that Eq. ( 12 ) has been 
obtained with the help of the relation 



At 



m 



m 



t-2 



2m, 



sin(20 t ) + /i susy cot (3 . 



(13) 



where A t is the soft SUSY breaking trilinear coupling of the top squarks. 

In the numerical results which are presented in Section [3] we include also effects from a 
non-vanishing bottom quark mass. The corresponding formulae are taken over from the 



literature 19,20 and are included in our program which computes the cross section. Thus, 
up to NLO we can study the production of the light MSSM Higgs boson also for relatively 
small values of Ma and large tan /3, the region where the bottom quark contributions are 
numerically important. At NNLO only those corrections are included which survive in 
the limit m& — > 0. 



2.3 Constructing approximation results from expansions 

We evaluate the three-loop corrections in the following hierarchies 

(hi) mg w mfj ~ m£ 2 ~ m-g 3> m t , 



11 



(h2) rrig pa m t ~ 2 pa m g 3> mi 3> m t , 

(h3) rriq pa mj 2 ~ m g ^> mi x ~ m t , (14) 

where m 4 , and represent the top quark, top squark, and gluino mass, respectively. 
trig is the generic mass of the supersymmetric partners of one of the five light quarks which 
we assume degenerate. Thus we only allow for a mixing in the top squark sector. A strong 



hierarchy between masses in Eq. ( 14 ) suggests that we apply an asymptotic expansion (see 



case an approximation sign is present. At three-loop order terms up to 0(l/m^ nsy ) have 
been computed for (hi) and (h3)j^ and up to 0(l/m^ ) and 0(l/m^ usy ) for (h2). For each 
mass difference at least three expansion terms (i.e. terms including (m 2 — m 2 ) 2 ) have 
been evaluated. 

In this context let us mention that in the hierarchies of Eq. ( JT4| ) we have omitted the e 
scalar mass which is simultaneously integrated out in order to guarantee the immediate 
matching to the Lagrangian in Eq. ([!]) regularized using DREGj^] We assume M £ to 
be much smaller than all other heavy masses which facilitates the calculation. It is 
furthermore consistent with the condition that the on-shell renormalized e scalar mass 
vanishes (see above). 

Experience shows that the Taylor expansion around equal masses works quite well and 
a few terms are sufficient in order to have good approximations for \m\jmi — 1| < 0.5 



(see, e.g., Refs. 31 43] for applications within a similar framework). This is the main 
motivation for hierarchy (hi) which covers a big part of all interesting cases. We have 
computed C\ for the hierarchy (h2) in order to be able to cover the case where one of the 
top squarks are light whereas the remaining supersymmetric masses are above 1 TeV and 
close together. The motivation for (h3) comes from situations of (h2) where m^ and m t 
have the same order of magnitude which is taken care of by a Taylor expansion in the 
mass difference. 



The strong hierarchies in Eqs. (14) do not leave any freedom in view of the expansion 
parameter which is essentially given by the inverse heavy mass. The choice of the mass 
in the numerator is practically irrelevant. This is different for the expansions in mass 
differences. There is a vast variety of choices which can lead to significant differences 
in the final result. In the following we classify the various possibilities and construct a 
procedure to obtain an optimized prediction. 



The results for the hierarchies of Eq. ( 14 ) contain either linear or quadratic differences of 



the involved masses, mi, mi ,m^ and rrig, which we define as 

= mi- mj , 

Ag = m 2 -m 2 , (15) 
respectively, and we consider expansions in A^/m, and Afj/m 2 . For the evaluation of the 



1 m susy is the common heavy mass in these hierarchie s 



In this context we also refer to the works 29 40 42 
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Figure 4: Figures (a) to (c) show different ways of expanding in small mass differences, 
when the masses rrij, rrij and m k are of the same order of magnitude. In (a) the expansion 
is performed around the basis mass mj, in (b) around rrij and in (c) around m k . Figure (d) 
demonstrates the expansion around a generic reference mass rriR. For the special values 
of rriR equal to mj,mj, or m k the representations of (a), (b) and (c) are reproduced. 



integrals in the individual hierarchies we have chosen the following differences: 

0*)'- ^.A^A*. (16) 

Note, however, this particular choice of representation may not be well suited for the 
numerical evaluation of certain mass spectra. For example, if in the case of (hi) m~ tx 
is the lightest supersymmetric mass one observes a worse convergence as compared to 
a representation where the expansion is carried out around one of the heavier masses, 
because higher order terms in the original expansion are less suppressed. 

The choice of the expansion point for the mass differences is akin to the choice of the 
renormalization scale /xr. Often it is argued that good choices go along with small higher- 
order corrections since certain classes of corrections are automatically summed. One could 
imagine a similar mechanism when expanding in mass differences: a clever choice of the 
basis mass can lead to small higher order terms in the expansion parameter. This criterion 
will be used below for finding the best suited representation. 

We have worked out two different approaches to optimize the proper choice of the expan- 
sion parameters in an automatic way. 

In the first approach all possible representations of the underlying expansion are generated 
automatically both for linear or quadratic mass differences A^ and A^. To this end in 
a first step the existing representation is rewritten in terms of new mass differences for 
each possible choice of the basis mass which is illustrated in Fig. |4](a)-(c). Afterwards an 
expansion in small quantities is performed respecting the original expansion depths. This 
approach has already been introduced in Ref. (44]. 

For the hierarchy (hi) we have four different basis masses {mg,mi ,mi 2 ,m%) and two 
possibilities (linear, quadratic) to write each of the three mass differences. Thus one has 
4 x 2 3 = 32 different representations. For the hierarchy (h2) we have 3 x 2 2 = 12 different 
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representations and for (h3) one obtains (2 x 2 1 ) x (3 x 2 2 ) = 48 different representations 
where the first bracket results from the expansion m t « m tl and the second one from 
m- t2 rs m,j « rrig. 

After generating all representations for a certain hierarchy they are evaluated numerically 
for a given supersymmetric spectrum. In order to select the representation which leads 
to the best approximation we define two quantities which we require to be small (see also 
Ref. |29|). The first one controls the deviation to the exact two-loop result and is defined 
through 

c (l),app _ (l),exact 

61 = (ij^tt ' ( 17 ) 

1 



c 



where c^' app contains as many expansion terms as the three-loop result we are interested 
in. The second quantity is given by 



(2),app,c _ (2), app 

Oo = — 



~ C l 

S 2 )' a pp 
-1 



where the superscript "c" indicates that the expansion is mass differences is reduced from 
three to two terms. 62 is a measure for the importance of higher order terms and thus the 
convergence of the expansions. 

The best representation is chosen by minimizing the combination 

5app= VW+W, (19) 

i.e., we require both good agreement with the exact result at two loops and a good 
convergence within the three-loop expression. 

Our second method is a generalization of the first one in the sense that for the basis 
mass we do not use one of the physical masses involved in the calculation but a generic 
reference mass tur. Of course, the numerical value of rriR should be of the order of the 
other masses, however, apart from that there is no restriction. This situation is illustrated 
in Fig. |4^d). Choosing a variable (although definite) basis mass reduces the number of 
representations, Nr, which can be seen from the following formulae 

2 max 

N ( R m) = U" 2 
i=i 

^ max 

4 bm) = n 2 " ( 2 °) 
1=1 

where "fbm" (fixed-basis-mass) and "vbm" (variable-basis-mass) refers to the first and 



second approach, respectively. In Eq. (20) i max refers to the number of independent 



regions with expansions in mass differences (for (hi) and (h2) we have i max = 1 and for 
(h3) z max = 2) and rii to the number of physical particle masses involved in the i-th region. 
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In the "vbm" scheme this leads to 16, 8 and 32 different representations for the hierarchies 
(hi), (h2) and (h3), respectively. 

Let us next discuss the criteria which fix m R in the "vbm" scheme. Note that the exact 
two- loop result is independent of m R and thus m R can be chosen in such a way that 
the approximated two- loop expression has minimal deviation from the exact result. This 



means that one has to find the minimum of the function |5i(m^)| (see Eq. (17)) which 



leads to the reference mass m R . The latter can be used to obtain the three-loop term 



of the matching coefficient, c\ 



R 

( 2 V™(ft) 



m 



R 



Altogether we have 



(ft) 



\5i(m R )\=mm => m R ^ V ,„ R 

where for the starting value of m R we choose the average over all involved masses. 



In complete analogy we determine the quantities m R 



(-52 ) 



m 



(5 1+2 ) (2), (5 2 ) 



(2), (5i +2 ) 



R 



1 



R 



(21) 



and 



i 



R 



Y 



\^2{m R )\ 
\Si(m R )\ + \5 2 (m R )\ 



mm 



mm 



rn 



m 



(ft) 

R 

(Sl+2) 

R 



C 



(%™(ft) 
R 



III 



C 



(2)/ (ft+2). 



[m 



R 



(22) 



For well-converging expansions we expect similar values of m R for all three choices of 5. 
The information collected so far is used to compute the uncertainty of the approximation 

(2) 

for c\ . We define it to be proportional to the deviations of the approximated results 



based on the criteria in Eqs. (21) and (22) to their average. Since this term could by 



chance be small an additional factor is introduced which takes into account the size of 5\ 
and 62- Our final formula for the uncertainty reads 



(2) 



1 + 



where is given by: 



,(2) 



+ 



<5=<5i ,52,(5i+2 



<5=<5i,<52,<5i+2 



(2) A™ W 
R 



III 



,(2) 



(23) 



(24) 



The described procedure leads to predictions for the three-loop coefficient in the form 



„(2), 



± Ac\ (i) where i runs over all representations. For our final result we pick the one 
with the smallest error, after ensuring that all predictions agree within their errors. In 
order to fulfill the latter step it may happen that the uncertainties have to be rescaled. 



Note that in all results presented in Section |3| Acjj 2 ^ is negligibly small. 



2.4 C\ to two loops 

In this subsection we discuss the application of the formalism described above to the 
two-loop result for C\ where the occurring momentum integrals have been evaluated both 
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Figure 5: 

matrix as a function of the soft SUSY breaking parameter rh t . 



mi and 9t obtained from the diagonalization of the top squark mass 



by taking into account the exact dependence on all mass parameters and by applying the 



hierarchies defined in Eq. (14). A detailed numerical discussion of the cross section will 
be presented in Section |3j In the following we want to test the approximation procedure 
and discuss its accuracy. For this purpose we adopt the following values for the input 
parameters 

m p = 155 GeV , A t = 100 GeV , m Qa = 500 GeV , 

/i susy = 100 GeV , tan /3 = 10 . (25) 

where A t is the trilinear coupling, Mz is the Z boson mass and rfiQ 3 , a soft SUSY breaking 
parameter for the squark doublet of the third family. We furthermore fix the renormal- 
ization scale to = M t , where M t is the on-shell top quark mass. 

Once these parameters are fixed it is possible to compute the physical masses , m ( - 2 and 
9t in dependence on the parameter rht, the singlet soft SUSY breaking mass parameter of 
the right-handed top squarkj^] This is achieved by diagonalizing the corresponding mass 
matrix. For convenience we show in Fig. [5] , m^ 2 and 8 t for m t between and 1 TeV. 
For small and large values of fh t one obtains a large splitting in the top squark masses 
whereas for m t f» 500 GeV they are approximately degenerate. 

The numerical result for is shown in Fig. ^ which is obtained by using m- t2 = nig = 
rrig and a = —0.1. Note that the hierarchies (h3) and (h2) are designed to lead to 
good approximations for small and large rht whereas (hi) should lead to good results for 



See Ref. 45 for details on the notation. 
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Figure 6: Two-loop correction to C% as a function of fh t . See text for further details. 



intermediate values. In Fig. [6]the exact result is shown as solid line and the approximations 
for the hierarchies (hi), (h2) and (h3) as dashed lines, where we also plot the corresponding 



uncertainties obtained from the "vbm" method described in Subsection 2.3 For clarity. 



the curves are truncated beyond the points where the corresponding uncertainty grows 
too large. The dotted lines correspond to the "fbm" method. They are different from the 
dashed ones only for a few values of fh t . From Fig. [6] one observes that for each value 
of fht there is at least one hierarchy which leads to an excellent approximation of the 
exact curve. In particular, good agreement is observed in the overlapping regions around 
fn t pa 300 GeV and 600 GeV. Furthermore, the hierarchy (h2) leads to good results also 
for small values of fn t which is not expected in the first place. It seems, however, that 
the expansion for 3> m t converges well. Moreover, one can learn from Fig. [6] that the 
approximation procedure seems to work very well since the estimated uncertainties cover 
the differences between the different hierarchies. 

Let us briefly comment on the limit m t — > where q 1 becomes numerically large. It 
originates from contributions containing self energy corrections to the t\ where in the loop 
heavy supersymmetric masses are present. In the DR scheme this leads to enhancement 
factors proportional to m^ usy / which become large for small values of . As a way out 



it is possible to use the MDR scheme introduced in Ref. 31 which systematically subtracts 



these terms. In the phenomenological results presented in Section [3] the supersymmetric 

^susy / 



parameters are such that the factors /m~ are not numerically dominant and thus 



we are not forced to introduce the MDR scheme. 
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Figure 7: Three- loop correction to C\ as a function of rh t . See text for further details. 



At this point a detailed comparison to evalcsusy is in orderf] In Ref. (23] DRED has 
been implemented in a naive way by performing the Lorentz algebra in four and the loop 
integrals in D = 4 — 2e dimensions, no e scalars have been introduced. Furthermore, 
evalcsusy only contains contributions where the Higgs boson couples to a top quark or 
top squark, diagrams as the ones shown in Fig. [2] are not included. Thus the result for 
C\ contains terms of order 1, m^/m-, ii susy m t /m'~ and Mf/m~ (as can also be deduced 
from the structure of the couplings as displayed in the Appendix of Ref. [23]). As already 
mentioned above, the M§/m| term is numerically small. The NLO calculation performed 
in this paper reproduces all terms contained in evalcsusy. 

Let us stress that it is not possible to apply the naive approach of Ref. [23] at NNLO. 



2.5 C\ to three loops 



In this subsection we use the parameters of Eq. (25) and the procedure described below 
this equation to evaluate the three- loop corrections to G\. The result is shown in Fig. [7] 
where the same notation as at two loops (see Fig. |6| has been used. The qualitative 



behaviour of the three-loop curves is very similar to the corresponding two-loop results 
in Fig. [6] Again one observes a nice overlap between the different hierarchies such that 
the whole parameter region can be covered. 

The analytical expressions are again quite lengthy and thus we refrain from presenting 
all of them in the paper but provide computer-readable expressions 38 



To illustrate the 



4 Note that the authors of Ref. 24 found agreement with the results of Ref. 23 . 
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structure of the results we list C\ for the hierarchy (hi) in the degenerate mass limit. It 
reads 



DR,(hl) 
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(26) 



As compared to Ref. [3] the Mf terms, which result both from the top quark/top squark 
and light quark/squark sector, are new. 



3 Numerical predictions 
3.1 Formalism 

As far as the contribution to the cross section originating from the top quark, top squark 
and squark sector is concerned we adopt the formalism outlined in Ref. [3] . In particular, 
we apply a redefinition of the effective-theory operator in such a way that both the effective 
operator and the redefined Wilson coefficient are renormalization group invariant. This 
allows us to separate the renormalization scale into a hard part, which enters the matching 
coefficient, and a soft part present in the matrix elements computed with the effective 
operator. In this way potentially large logarithms ln(/x^/M/j) are resummed where /i^ is of 
the order of the top quark or SUSY particle masses. In this framework the cross section 
can be computed from the following formula: 

af QCD (/^) = <4s (0) + ^-^E« + ^^2c«£(°> + ( ) E< 2 ) (27. 

7T 7T y \ 71 
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+ 



a^jfis) ^ 5 Vfe) 2c g) £ (i) + f as 5) (Vh) \ j( c (i))2 + 2c ( 2 )j E (o) 



+ . . . 



where 



l + «gVO c(1) + / \ c ( 2) 



(5), 



c 



+ 



7T 



7T 



E = E (o) , «gVg) s (i) , ( c*s 0, M \ S ( 2) 



(5), 



+ 



7T 



7T 



(28) 



and 



(TO 



: -I(r t ) +^2Si< —c et s et ( t a + 



2 2 

2 *i t2 AQ 

b 26 t 2 



+ s 



"+/3q 2~ 



c 26»t ( - 1 ] + S; 



M 2 ~ 



(29) 



with 



Am 2 



, Si — — 1 , S*2 — +1 , 



J(r) =r[l + (1 - r)/(r)] , J(r) = r(l - r/(r)) , 
/(r) = 



-4 In 



arcsin 2 (l/-y/r), r > 1 
Z7T ) , r < 



1 (1„ l + Vl-r 



(30) 



The coefficient function C g is obtained from the renormalized Wilson coefficient C\ in 
Eq. Q via 



G, 



3tt 1 

joy 5(6) 



Ci 



(31) 



with 



S (5) 



7T 



!/3(5) 



(5)^,(5) 



(32) 



Eq. (27 



In E q. (|31[ ) we divide by cf^ since it is already included in a which is factored out from 
Eq. (27). Note that the coefficients Cg and c g 2 ^ depend on fj, h which is suppressed in 



£( fc ) depends on /i s and represents the N fc LO contribution of the hadronic cross section 
computed within five-flavour QCD using the effective Lagrange density in Eq. where a 
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summation over all sub channels (gg, qg, qq, qq and qq ! ) is understood. Note that this part 
can be taken over from the corresponding SM calculation; in our case we use the results 
from Ref. 18 . In principle T,^ also depends on the factorization scale however, in our 
calculation we identify /if and /i s . £ is by construction renormalization group invariant 
since is incorporated. Note that for fXh = ^ s B^ cancels in the above formulae. 

The approach described so far is not practical in case the bottom quark mass is kept 
different from zero since then a division into fih and fi s is not as obvious as in the top and 
squark sector where the calculation is based on the effective theory. Thus the numerical 
evaluation of the cross section including finite bottom quark mass effects proceeds as 
follows: we set \ih = fi s and evaluate the cross section up to NLO using the expressions 
from Ref. Il9). Afterwards we subtract our NLO top- and squark-sector contribution 
based on Eq. (27) choosing also fih = fi s and finally add our NNLO result as given in 



Eq. (27). At the end we include electroweak corrections assuming complete factorization. 



Thus our best prediction looks as follows 

a SQCD {pp ^ h + x) = ( 1 + (J EW) x 



(33) 



a 



SQCD 
tb 



a 



SQCD 
t 



NLO 



+ <7. 



SQCD 



NLO 



NNLO 



where 5 EW is taken from Ref. [46] (see also Ref. 47 ). Note that in order to obtain 
the NNLO prediction the convolution of all contributions entering Eq. (33) is performed 
with NNLO parton distribution functions (PDFs). This is also true for the quantities 
°f& QCD (/-OInlo an d crf^ CD (Ais)|NLO- Analogous formulae to Eq. (33) are also used for LO 
and NLO predictions of a(pp — >■ h + X). Of course, in these cases LO and NLO PDFs 
are used, respectively. 



3.2 Numerical setup 

In the following we briefly describe the important features of our numerical program 
gghXsection. It consists of several components which take care of the generation of 
the spectrum, the computation of the Higgs boson mass, the convolution of the partonic 
cross sections with the PDFs in the effective theory and the determination of the best 
approximation for the three-loop matching coefficient. A flow chart showing the individual 
components can be found in Fig. |HJ 

• The individual components are called and controlled by the Mathematica program 
gghXsection. 

• Our analytical expressions are parametrized in terms of DR mass parameters de- 
fined at the renormalization scale /zr. For an automated generation of the required 
MSSM parameter spectra within Mathematica we use our own Mathematica pack- 
age LHSQLDBj^This package works as interface between Mathematica and the SUSY 

5 LHSQLDB stands for Les Houches Spectra & SQL Data Base, where the two S are melted to one. 
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Figure 8: Flow chart of our numerical setup. The rectangular frames denote stand-alone 
program packages which are combined by gghXsection. It controls the flow of data 
(indicated by the ovals) among the individual components. More details are given in the 
text. 



spectrum generators following the "Les Houches Accord' 

[52 



Spheno 51 



48,49 likeSOFTSUSY 50 



or Suspect |52| (we use SOFTSUSY as default spectrum generator 
throughout this paper). Furthermore it is able to connect Mathematica to an SQL 
data base where SUSY spectra can be automatically saved to and/or loaded from, 
avoiding multiple recalculations of already known parameters. 

Since this package might be of general use, a separate publication is planned. How- 
ever, a private version can be obtained from the authors upon request. 

We define the strong coupling constant in five-flavour QCD using the MS scheme. 
The numerical results for ai 5 \[iR) are taken from the PDF sets. 



The spectrum generator also provides input for H3m 31 53 which is used to compute 



the light MSSM Higgs boson mass to three-loop accuracy. H3m uses FeynHiggs [54] 
for the two-loop expressions. 

Note that the spectrum generators can in general only be used for renormalization 
scales as low as the Z boson mass. However, in our numerical analyses also lower 
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scales appear since we choose fi s = Mh/2 as a central value. For this reason we have 
implemented the renormalization group equations (taken from Ref. [55]) ourselves in 
order to obtain the spectrum for all desired values of n s (see the program RunMSUSY 
in Fig. D). 

In order to compute the DR top quark mass from the on-shell one we use the results 



of Ref. 56 ; 57 which are based on the program TSIL. For this step we need the strong 



coupling constant in the full theory which we obtain with the help of RunDec 58,59 



and the extension of the decoupling and running routines to supersymmetry based 



on Refs. 40,60 



The integrations needed for the computation of £ are performed with the help of 
the self- written C++ program Xsection. 



The SM part entering our numerical analysis is based on results obtained in Ref. 18 



It is used to extract the quantities £W in Eq. (27) which are necessary to construct 
the top sector result up to NNLO. 

In the C++ program we have also included the computation of cr^ QCD (/i s )| NL o which 
is based on the formulae provided in Refs. 19,23 . The analytic expressions for the 



amplitudes which we implemented in our numerical program are taken from Ref. [19] 
for the bottom and from Ref. 23 for the top sector which provides the SUSY QCD 



corrections up to NLO. Thus the exact bottom quark mass dependence is taken 
into account as well as the top quark, gluino and squark dependence in the effective- 
theory approach, i.e. in the formal limit Mh — > 0. For Higgs boson masses considered 
here this approximation approaches the exact result to high accuracy 15 18 . We 

in 
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have cross checked our implementation against the results implemented in 
certain limits. 



We use LO, NLO and NNLO PDFs by the MSTW2008 collaboration [6l] in order 



to obtain the respective predictions for cr s ^ CD in Eq. (33). The choice of the PDFs 
determines the values of a s (M z ) = ai 5 ^(M^): 



0.1394, c^ LO (M z ) = 0.1202, 



a 



NNLO 



{M z ) = 0.1171. (34) 



The appropriate MS beta function then determines a s (nn) that enters the formulae. 
The procedure for the calculation of the best three-loop approximation of C\ within 



SQCD as described in Subsection 2J3 is implemented in Mathematica. 

The electroweak correction which enter the final result in factorized form (cf. 



Eq. (33)) are taken from Ref. 46 



If not stated otherwise the following default values are used: 

^ = ^, »h = M t , M t = 173.3GeV. 



(35) 
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3.3 Prediction to NNLO 



In order to discuss the numerical effects we start with the m™ ax 



scenario introduced in 



Ref. 62 . The low-energy parameters are given by 

A'-susy 



A t = A b = A T = V6m susy + , 

tclR p 

M 1 = 5s 2 w /(3c 2 w )M 2 . 



A'-susy 

M 2 



200 GeV . 
200 GeV . 



m 



susy 



1000 GeV 



800 GeV , (36) 



where the mass of the pseudo-scalar Higgs boson Ma and the ratio of the vacuum ex- 



pectation values of the two Higgs doublets tan/3 is varied. In Eq. (36) A t is the trilinear 
coupling, /Xsusy is the Higgs-Higgsino bilinear coupling from the super potential, m 



"susy 



IS 



the common mass value for all squarks, Mi, M 2 and M 3 are the gaugino mass parameters 
and % and c\y are the sine and cosine of the weak mixing angle. The parameters in 



Eq. (36) are defined at the scale //r 



///-„.,„ + M t 2 . 



7,2 

susy 



In Fig. [9] we discuss the dependence of the cross section on tan/3 and Ma, the remaining 



parameters are fixed to the values of Eq. (36). The NNLO prediction for the top sector 
contribution cr t S< ^ CD (/i s , /i/j) is shown in Figf^a) where the numbers on the contour line^] 

6 The small gray numbers in the background provide the cross section in pb for various values of tan f3 
and Ma- In order not to overload the plot it was chosen to keep these numbers small which means 
that they probably can only be read after magnifying the plots on the screen. Nevertheless we prefer to 
provide this information. 
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Figure 10: (a) NNLO K factor corresponding to the prediction in Fig. [9]^b). Panel (b) 
shows the ratio of the MSSM and the SM NNLO prediction from the top quark/top squark 
sector. 



correspond to the cross section in pbQ Ma is varied between 100 GeV and 200 GeV and we 
choose tan (5 < 20. The (red) hyperbola-like lines indicate Higgs boson masses of 123 GeV 
(lower) and 129 GeV, respectively. Thus the area in between is in accordance with the 
result reported by ATLAS [I] and CMS The NNLO cross section has been obtained 
with the help of hierarchy (hi) which is appropriate since the SUSY mass spectrum is 
almost degenerate. 

For comparison we show in Fig. [9Tb) the complete result including also bottom quark 



effects and the electroweak corrections according to Eq. (33). One observes a significant 
enhancement of the cross section for small values of Ma and large tan /3 (i.e. in the 
upper left corner). In the remaining parts of the plot only a reduction of a few percent is 
observed which comes basically from the destructive interference of the top and bottom 
contribution. Let us mention that a large part of the parameter space shown in Fig. [9] is 
already excluded experimentally which will be discussed below. Basically only the lower 
right part of the plot is experimentally still allowed. Note that this region is dominated 
by the contribution from the top sector and the bottom contribution only plays a minor 
role. 

The effect of higher order corrections is conveniently parametrized in terms of the K 



7 Note that er^ CD contains also contributions where the Higgs boson couples to squarks corresponding 
to light quarks. Since these contributions are numerically small we write in the legend only u t + t" . 
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factor which we define through 



a HO 



K = ^LO- ( 37 ) 



In Fig. [nja) we show the NNLO (i.e. HO=NNLO in Eq. ([37])) K factor for the same 
parameter space as in Fig. [9] In the experimentally allowed region it amounts to about 
100% where about 70% arise from the NLO corrections. It is interesting to compare the 
MSSM prediction to the one of the SM which is done in Fig. [To](b) where the ratio of the 
NNLO predictions in the two models is plotted. One observes that the supersymmetric 
corrections reduce the cross section by typically 60% to 90% in the parameter space 
where the Higgs boson mass is between 123 GeV and 129 GeV. Note, however, that for 
Ma ^ 140 GeV this effect basically results from a modification of the Higgs couplings to 
fermions and thus is already present at LO. In fact for this region the K factors in the 



SM and MSSM are very similar for the m™ ax scenario. In Fig. 10 b) we also show the 
regions excluded by searches performed at LEP [63) and CMS [64] which is indicated by the 
gray areaj^] 

The reason for the small deviations between the SM and MSSM in the m™' 



the heavy supersymmetric mass spectrum resulting from the input parameters in Eq. (36). 
Actually, the typical squark mass is 1000 GeV with only small differences between the 
individual flavours and the gluino mass takes the value m- g = 860 GeV. 

In the following we slightly modify the parameters in such a way that one of the top 
squarks becomes light whereas the other remains heavy. Thus, on one hand one obtains 
sizeable supersymmetric contributions to the Higgs boson production cross section and 
at the same time a Higgs boson mass which is sufficiently heavy. To date such a scenario 
has not been excluded by experimental searches. 

In our modified m™ ax scenario the default values of the low-energy parameters are given 
by 

A b = A T = 2469.48 GeV , A t = 1500 GeV , /i susy = 200 GeV , 
Mi = 5s 2 w /(3c 2 w )M 2 , M 2 = 200 GeV , M 3 = 800 GeV , 

m susy = 1000 GeV , m t = 400 GeV , 

M A = 1000 GeV , tan /3 = 20 , (38) 

where m t represents the singlet soft SUSY breaking parameter of the right-handed top 



squark (i.e. EXTPAR 46 in the "Les Houches Accord" notation 48,49]). 

In Fig. [IT[a) we show the top squark masses as a function of fht- The variation of this 
parameter essentially changes mi whereas mr remains almost constant at a large value 
of about 1050 GeV. For large values of rh t one obtains a heavy spectrum and thus the 
supersymmetric corrections to the Higgs boson production become small. On the other 
hand, for a small value of m t and thus a large splitting in the top squark sector the 



The data points have been extracted from Fig. 3 of Ref. 64 using the program EasyNData 65 
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Figure 11: (a) Top squark masses and (b) lightest Higgs boson mass as a function of fh t 



for the parameters given in Eq. (36). 



numerical effect can be large as we will discuss below. The hierarchy in the top squark 
sector suggests to use (h2) for evaluating the three-loop corrections to C\. 

It is interesting to have a closer look at the numerical values of the three-loop coefficient 
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cf } . For rht = 400 GeV we hav|] cj 2) (/i = M t ) « -0.58(1) and cf } (^ = M h /2 « 
63 GeV) w -14.66(1) and for ra t = 600 GeV one obtains cfV = M t ) « 5.70(1) and 
c[\n = M h /2 ?a 63 GeV) ~ —1.29(23). One observes a quite strong dependence of 

(2) 

c\ on the renormalization scale which demonstrates the importance of the three-loop 
calculation. Note that the renormalization scale dependence is compensated by the other 



quantities entering formula (27) for the cross section. 

In order to get a feeling about the physical masses we list in the following the DR mass 
values for the squark and gluino masses. For = M t and fh t = 400 GeV one obtains 
from the parameters provided by S0FTSUSY the following values 

m k = 370 GeV , m h = 1045 GeV , 

m q = 1042 GeV , m- g = 860 GeV , (39) 

where trig corresponds to the average of mj, toj, m s and mg. 



The lightest MSSM Higgs boson mass Mh as predicted from the parameters in Eq. (36) 



is shown in Fig. [TT](b) where m t is varied in the same range as before. In this figure we 



compare M h as computed by S0FTSUSY (dashed lines) to the one obtained by H3m 31 53 
(solid line). Note that the latter, which includes three-loop corrections, is more than 
1 GeV above the former. As one can see, in the whole range of fht the prediction remains 
within the horizontal lines which mark 123 GeV and 129 GeV, respectively, the limits 
motivated by the recent results reported by ATLAS [l] and CMS [2]. 



In Fig. 12 we show a t (fj, s , fih) as a function of fht at LO, NLO and NNLO (from bottom 
to top). For each order three curves are shown where the dotted curve corresponds to the 
SM. The SQCD corrections are included in the dashed and solid line where for the former 
the soft and hard scale have been identified with Mh/2 and for the latter fi s = Mh/2 
and fih = M t has been chosen. As already mentioned above, the difference between the 
SM and MSSM corrections is small for large values of fht, for smaller values, however, a 
sizeable effect is visible. For example, for fh t = 400 GeV a reduction of the cross section 
of about 5% is observed. 

The difference between the dashed and solid line quantifies the effect of the resummation 



discussed in Subsection 3.1 It is negligible for large fh t , however, for smaller values it can 
be as large as 2-3%. It is interesting to note that supersymmetric three- loop corrections 
to C\ computed in this paper play an important role. In case the three-loop coefficient is 
identified with the SM one a reduction of only 3% and not 5% is observed. 



In Fig. [13] we fix fht = 400 GeV and discuss the renormalization scale dependence of the 
LO, NLO and NNLO result by varying y, h = fi s between 25 GeV and about 400 GeV. The 
SM (dotted) and the MSSM (solid lines) curves show a very similar behaviour. Close to 
the central scale fi s = Mh/2 one observes a crossing of the NNLO and NLO curves, i.e., 
vanishing NNLO corrections. It is at a slightly higher value for the MSSM which can be 
explained by the larger masses in the spectrum entering the predictions. 

9 The numbers in the brackets indicate the uncertainty according to our approximation procedure. 
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Figure 12: The cross section <7 f SQCD (/i s , /i/j) as a function of m t using the parameters in 
Eq. (36) to LO (bottom), NLO (middle) and NNLO (top). The dotted line corresponds 



to the SM and the dashed and solid lines to the MSSM. 
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Figure 13: Renormalization scale dependence of the cross section <7 f SQCD (yU s 



fih)- 



In a next step we study the dependence on tan j3 for fixed m t = 400 GeV. One observes 
top squark masses which are almost constant and are given by m t - f» 370 GeV and 
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Figure 14: Lightest Higgs boson mass as a function of tan for the parameters given in 
Eq. (ph. 



mi 2 w 1040 GeV. On the other hand there is a strong dependence of Mh on tan j3 which 
is shown in Fig. 14 for the result from SOFTSUSY and H3m. As one can see, for tan /3 < 6 
the predicted Higgs boson mass is below 123 GeV using H3m; on the basis of SOFTSUSY 
this boundary is reached already for tan /3 ~ 9. 



In Fig. 15 we show o"f^ CD (/i s , fih) as a function of tan/3 for 5 < tan/3 < 30 at LO, NLO 
and NNLO (from bottom to top) choosing A t = 1500 GeV and rh t = 400 GeV. The 
SM result is again shown as dotted line. The MSSM predictions are both shown for 
/J>h = /J-s = Mh/2 (dashed) and \ih = M t and /x s = Mh/2 (solid line). As compared to 
the SM prediction the MSSM results are again reduced by a few percent. Note that this 
effect increases when going from LO to NLO and finally to NNLO where a difference of 
about 5% is observed. It is also interesting to mention that the proper choice of scales is 
not negligible: the results for fih = are in general a few percent above the ones with 
fJ-h ^ /V 

Let us finally present results for a(pp — » h + X) which include in addition bottom quark 
contributions up to NLO and furthermore also electroweak corrections. In Fig. [XB|a) 
and (b) the dependence on m t and tan /3 is plotted, respectively, where from bottom to 
top the LO, NLO and NNLO results are shown. We refrain from showing the bottom 
quark contributions separately since they are about two orders of magnitude smaller than 
the top contributions. However, the interference terms have a visible effect which can 
be seen by the difference between the dashed and dotted curves since in the latter the 
quantities <7^ CD (/i s )| NL o and cr^ CD (/i s )| NLO in Eq. (33) are set to zero. One observes a 
reduction of about 5% after including bottom quark effects — basically independent of rh t 
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Figure 15: The cross section crf QCD (/i s , fih) as a function of tan/3 using the parameters in 



Eq. (38) to LO (bottom), NLO (middle) and NNLO (top). The dotted line corresponds 
to the SM and the dashed and solid lines to the MSSM. 



and tan f3. Thus even for tan f3 = 20 the bottom quark effects are small for the considered 
scenarios and hence a NNLO calculation for m;, ^ is not mandatory. The reduction from 
top/bottom interference is to a large extend compensated by the electroweak corrections 
taken into account in a multiplicative way (cf. Eq. (33)) as can be seen by solid line which 
includes all contributions of Eq. (33). 



4 Conclusions 

We have computed the three-loop corrections to the matching coefficient of the effective 
Higgs-gluon coupling C\ originating from supersymmetric QCD. The three-loop integrals 
contain several mass scales and thus an exact calculation on the basis of the techniques 
available these days is not possible. We have developed an approximation method based on 
expansions in mass hierarchies and mass differences which selects the best parametrization 
for the expansion parameters. The described procedure is not bound to the problem at 
hand but can certainly also be applied to other processes. 

The three-loop result for G\ is used in order to evaluate the Higgs boson production 
cross section at LHC to NNLO. Whereas the three-loop expressions are only available 
for vanishing light quark masses we include up to NLO the contributions from the top 
and bottom sector allowing for a non-vanishing bottom quark mass. Numerical results 
for the total production cross section are presented for a center-of-mass energy of y/s = 
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Figure 16: The cross section cr s ^ CD (j>p — > h + X) as a function of (a) m t and (b) tan/3 
using the parameters in Eq. ( 38 ) . 



8 TeV for two MSSM scenarios which predict a light supersymmetric Higgs boson mass in 
accordance with the recent results from ATLAS and CMS. We study both the dependence 
on the renormalization scale, tan (5 and on the m t , the singlet soft SUSY breaking mass 
parameter of the right-handed top squark. For the case where all supersymmetric masses 
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are above 1 TeV the supersymmetric corrections are small. However, in case one of the top 
squarks becomes of the order of a few hundred GeV the supersymmetric corrections to the 
matching coefficient C\ are not negligible any more and a sizeable effect on the production 
cross section is noticeable. For the scenarios considered in this paper a reduction of the 
NNLO prediction of the cross section by a few percent is observed. We also want to stress 
that visible effects result from separating the soft and hard scale in C\ which automatically 
resums potentially large logarithms involving the Higgs boson mass. 

In order to perform the numerical evaluation a flexible tool has been developed which 
allows a convenient handling of the spectrum generators in combination with Mathematica 
and C++ routines evaluating the cross sections. In a comfortable way it is possible to vary 
the renormalization scale or other parameters entering the prediction. 
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